Genotyping for genetic association studies: Methods and applications

In this thesis, two separate single nucleotide polymorphism (SNP) genotyping techniques were set up at the Finnish Genome Center, pooled genotyping was evaluated as a screening method for large-scale association studies, and finally, the former approaches were used to identify genetic factors predisposing to two distinct complex diseases by utilizing large epidemiological cohorts and also taking environmental factors into account. The first genotyping platform was based on traditional but improved restriction-fragment-length-polymorphism (RFLP) utilizing 384-microtiter well plates, multiplexing, small reaction volumes (5 µl), and automated genotype calling. We participated in the development of the second genotyping method, based on single nucleotide primer extension (SNuPeTM by Amersham Biosciences), by carrying out the alpha- and beta tests for the chemistry and the allele-calling software. Both techniques proved to be accurate, reliable, and suitable for projects with thousands of samples and tens of markers. Pooled genotyping (genotyping of pooled instead of individual DNA samples) was evaluated with Sequenom s MassArray MALDI-TOF, in addition to SNuPeTM and PCR-RFLP techniques. We used MassArray mainly as a point of comparison, because it is known to be well suited for pooled genotyping. All three methods were shown to be accurate, the standard deviations between measurements being 0.017 for the MassArray, 0.022 for the PCR-RFLP, and 0.026 for the SNuPeTM. The largest source of error in the process of pooled genotyping was shown to be the volumetric error, i.e., the preparation of pools. We also demonstrated that it would have been possible to narrow down the genetic locus underlying congenital chloride diarrhea (CLD), an autosomal recessive disorder, by using the pooling technique instead of genotyping individual samples. Although the approach seems to be well suited for traditional case-control studies, it is difficult to apply if any kind of stratification based on environmental factors is needed. Therefore we chose to continue with individual genotyping in the following association studies. Samples in the two separate large epidemiological cohorts were genotyped with the PCR-RFLP and SNuPeTM techniques. The first of these association studies concerned various pregnancy complications among 100,000 consecutive pregnancies in Finland, of which we genotyped 2292 patients and controls, in addition to a population sample of 644 blood donors, with 7 polymorphisms in the potentially thrombotic genes. In this thesis, the analysis of a sub-study of pregnancy-related venous thromboses was included. We showed that the impact of factor V Leiden polymorphism on pregnancy-related venous thrombosis, but not the other tested polymorphisms, was fairly large (odds ratio 11.6; 95% CI 3.6-33.6), and increased multiplicatively when combined with other risk factors such as obesity or advanced age. Owing to our study design, we were also able to estimate the risks at the population level. The second epidemiological cohort was the Helsinki Birth Cohort of men and women who were born during 1924-1933 in Helsinki. The aim was to identify genetic factors that might modify the well known link between small birth size and adult metabolic diseases, such as type 2 diabetes and impaired glucose tolerance. Among ~500 individuals with detailed birth measurements and current metabolic profile, we found that an insertion/deletion polymorphism of the angiotensin converting enzyme (ACE) gene was associated with the duration of gestation, and weight and length at birth. Interestingly, the ACE insertion allele was also associated with higher indices of insulin secretion (p=0.0004) in adult life, but only among individuals who were born small (those among the lowest third of birth weight). Likewise, low birth weight was associated with higher indices of insulin secretion (p=0.003), but only among carriers of the ACE insertion allele. The association with birth measurements was also found with a common haplotype of the glucocorticoid receptor (GR) gene. Furthermore, the association between short length at birth and adult impaired glucose tolerance was confined to carriers of this haplotype (p=0.007). These associations exemplify the interaction between environmental factors and genotype, which, possibly due to altered gene expression, predisposes to complex metabolic diseases. Indeed, we showed that the common GR gene haplotype associated with reduced mRNA expression in thymus of three individuals (p=0.0002).

Statistical analysis of the associations of hereditary factors with the risk of Type 1 diabetes and Diabetic Nephropathy based on the familial data collected through ascertaiment from population-based registers

In genetic epidemiology, population-based disease registries are commonly used to collect genotype or other risk factor information concerning affected subjects and their relatives. This work presents two new approaches for the statistical inference of ascertained data: a conditional and full likelihood approaches for the disease with variable age at onset phenotype using familial data obtained from population-based registry of incident cases. The aim is to obtain statistically reliable estimates of the general population parameters. The statistical analysis of familial data with variable age at onset becomes more complicated when some of the study subjects are non-susceptible, that is to say these subjects never get the disease. A statistical model for a variable age at onset with long-term survivors is proposed for studies of familial aggregation, using latent variable approach, as well as for prospective studies of genetic association studies with candidate genes. In addition, we explore the possibility of a genetic explanation of the observed increase in the incidence of Type 1 diabetes (T1D) in Finland in recent decades and the hypothesis of non-Mendelian transmission of T1D associated genes. Both classical and Bayesian statistical inference were used in the modelling and estimation. Despite the fact that this work contains five studies with different statistical models, they all concern data obtained from nationwide registries of T1D and genetics of T1D. In the analyses of T1D data, non-Mendelian transmission of T1D susceptibility alleles was not observed. In addition, non-Mendelian transmission of T1D susceptibility genes did not make a plausible explanation for the increase in T1D incidence in Finland. Instead, the Human Leucocyte Antigen associations with T1D were confirmed in the population-based analysis, which combines T1D registry information, reference sample of healthy subjects and birth cohort information of the Finnish population. Finally, a substantial familial variation in the susceptibility of T1D nephropathy was observed. The presented studies show the benefits of sophisticated statistical modelling to explore risk factors for complex diseases.

Genetic variation in the LDL receptor-related protein 5 (LRP5) gene : Association with bone health and metabolic parameters

Bone mass accrual and maintenance are regulated by a complex interplay between genetic and environmental factors. Recent studies have revealed an important role for the low-density lipoprotein receptor-related protein 5 (LRP5) in this process. The aim of this thesis study was to identify novel variants in the LRP5 gene and to further elucidate the association of LRP5 and its variants with various bone health related clinical characteristics. The results of our studies show that loss-of-function mutations in LRP5 cause severe osteoporosis not only in homozygous subjects but also in the carriers of these mutations, who have significantly reduced bone mineral density (BMD) and increased susceptibility to fractures. In addition, we demonstrated for the first time that a common polymorphic LRP5 variant (p.A1330V) was associated with reduced peak bone mass, an important determinant of BMD and osteoporosis in later life. The results from these two studies are concordant with results seen in other studies on LRP5 mutations and in association studies linking genetic variation in LRP5 with BMD and osteoporosis. Several rare LRP5 variants were identified in children with recurrent fractures. Sequencing and multiplex ligation-dependent probe amplification (MLPA) analyses revealed no disease-causing mutations or whole-exon deletions. Our findings from clinical assessments and family-based genotype-phenotype studies suggested that the rare LRP5 variants identified are not the definite cause of fractures in these children. Clinical assessments of our study subjects with LPR5 mutations revealed an unexpectedly high prevalence of impaired glucose tolerance and dyslipidaemia. Moreover, in subsequent studies we discovered that common polymorphic LRP5 variants are associated with unfavorable metabolic characteristics. Changes in lipid profile were already apparent in pre-pubertal children. These results, together with the findings from other studies, suggest an important role for LRP5 also in glucose and lipid metabolism. Our results underscore the important role of LRP5 not only in bone mass accrual and maintenance of skeletal health but also in glucose and lipid metabolism. The role of LRP5 in bone metabolism has long been studied, but further studies with larger study cohorts are still needed to evaluate the specific role of LRP5 variants as metabolic risk factors.

Examination of Genetic Components Affecting Human Obesity-Related Quantitative Traits

Obesity increases the risk for several conditions, including type 2 diabetes mellitus, cardiovascular disease, hypertension, osteoarthirits and certain types of cancer. Twin- and family studies have shown that there is a major genetic component in the determination of body mass. In recent years several technological and scientific advance have been made in obesity research. For instance, novel replicated loci have been revealed by a number of genome wide association studies. This thesis aimed to investigate the association of genetic factors and obesity-related quantitative traits. The first study investigated the role of the lactase gene in anthropometric traits. We genetically defined lactose persistence by genotyping 31 720 individuals of European descent. We found that lactase persistence was significantly correlated with weight and body mass index but not with height. In the second study we performed the largest whole genome linkage scan for body mass index to date. The sample consisted of 4401 twin families and 10 535 individuals from six European countries. We found supporting evidence for two loci (3q29 and 7q36). We observed that the heritability estimate increased substantially when additional family members were removed from the analyses, which suggests reduced environmental variance in the twin sample. In the third study we assessed metabonomic, transcriptomic and genomic variation in a Finnish population cohort of 518 individuals. We formed gene expression networks to portray pathways and showed that a set of highly correlated genes of an inflammatory pathway associated with 80 serum metabolites (of 134 quantified measures). Strong association was found, for example, with several lipoprotein subclasses. We inferred causality by using genetic variation as anchors. The expression of the network genes was found to be dependent on the circulatory metabolite concentrations.

Genetic analysis of chromosomal regions 2q33, 7q32 and 19q13 in multiple sclerosis susceptibility

Multiple sclerosis (MS) is an immune-mediated demyelinating disorder of the central nervous system (CNS) affecting 0.1-0.2% of Northern European descent population. MS is considered to be a multifactorial disease, both environment and genetics play a role in its pathogenesis. Despite several decades of intense research, the etiological and pathogenic mechanisms underlying MS remain still largely unknown and no curative treatment exists. The genetic architecture underlying MS is complex with multiple genes involved. The strongest and the best characterized predisposing genetic factors for MS are located, as in other immune-mediated diseases, in the major histocompatibility complex (MHC) on chromosome 6. In humans MHC is called human leukocyte antigen (HLA). Alleles of the HLA locus have been found to associate strongly with MS and remained for many years the only consistently replicable genetic associations. However, recently other genes located outside the MHC region have been proposed as strong candidates for susceptibility to MS in several studies. In this thesis a new genetic locus located on chromosome 7q32, interferon regulatory factor 5 (IRF5), was identified in the susceptibility to MS. In particular, we found that common variation of the gene was associated with the disease in three different populations, Spanish, Swedish and Finnish. We also suggested a possible functional role for one of the risk alleles with impact on the expression of the IRF5 locus. Previous studies have pointed out a possible role played by chromosome 2q33 in the susceptibility to MS and other autoimmune disorders. The work described here also investigated the involvement of this chromosomal region in MS predisposition. After the detection of genetic association with 2q33 (article-1), we extended our analysis through fine-scale single nucleotide polymorphism (SNP) mapping to define further the contribution of this genomic area to disease pathogenesis (article-4). We found a trend (p=0.04) for association to MS with an intronic SNP located in the inducible T-cell co-stimulator (ICOS) gene, an important player in the co-stimulatory pathway of the immune system. Expression analysis of ICOS revealed a novel, previously uncharacterized, alternatively spliced isoform, lacking the extracellular domain that is needed for ligand binding. The stability of the newly-identified transcript variant and its subcellular localization were analyzed. These studies indicated that the novel isoform is stable and shows different subcellular localization as compared to full-length ICOS. The novel isoform might have a regulatory function, but further studies are required to elucidate its function. Chromosome 19q13 has been previously suggested as one of the genomic areas involved in MS predisposition. In several populations, suggestive linkage signals between MS predisposition and 19q13 have been obtained. Here, we analysed the role of allelic variation in 19q13 by family based association analysis in 782 MS families collected from Finland. In this dataset, we were not able to detect any statistically significant associations, although several previously suggested markers were included to the analysis. Replication of the previous findings on the basis of linkage disequilibrium between marker allele and disease/risk allele appears notoriously difficult because of limitations such as allelic heterogeneity. Re-sequencing based approaches may be required for elucidating the role of chromosome 19q13 with MS. This thesis has resulted in the identification of a new MS susceptibility locus (IRF5) previously associated with other inflammatory or autoimmune disorders, such as SLE. IRF5 is one of the mediators of interferons biological function. In addition to providing new insight in the possible pathogenetic pathway of the disease, this finding suggests that there might be common mechanisms between different immune-mediated disorders. Furthermore the work presented here has uncovered a novel isoform of ICOS, which may play a role in regulatory mechanisms of ICOS, an important mediator of lymphocyte activation. Further work is required to uncover its functions and possible involvement of the ICOS locus in MS susceptibility.

Genetic background of bipolar disorder and related cognitive impairments

Bipolar disorder (BP) is a complex psychiatric disorder characterized by episodes of mania and depression. BP affects approximately 1% of the world’s population and shows no difference in lifetime prevalence between males and females. BP arises from complex interactions among genetic, developmental and environmental factors, and it is likely that several predisposing genes are involved in BP. The genetic background of BP is still poorly understood, although intensive and long-lasting research has identified several chromosomal regions and genes involved in susceptibility to BP. This thesis work aims to identify the genetic variants that influence bipolar disorder in the Finnish population by candidate gene and genome-wide linkage analyses in families with many BP cases. In addition to diagnosis-based phenotypes, neuropsychological traits that can be seen as potential endophenotypes or intermediate traits for BP were analyzed. In the first part of the thesis, we examined the role of the allelic variants of the TSNAX/DISC1 gene cluster to psychotic and bipolar spectrum disorders and found association of distinct allelic haplotypes with these two groups of disorders. The haplotype at the 5’ end of the Disrupted-in-Schizophrenia-1 gene (DISC1) was over-transmitted to males with psychotic disorder (p = 0.008; for an extended haplotype p = 0.0007 with both genders), whereas haplotypes at the 3’ end of DISC1 associated with bipolar spectrum disorder (p = 0.0002; for an extended haplotype p = 0.0001). The variants of these haplotypes also showed association with different cognitive traits. The haplotypes at the 5’ end associated with perseverations and auditory attention, while the variants at the 3’ end associated with several cognitive traits including verbal fluency and psychomotor processing speed. Second, in our complete set of BP families with 723 individuals we studied six functional candidate genes from three distinct signalling systems: serotonin-related genes (SLC6A4 and TPH2), BDNF -related genes (BDNF, CREB1 and NTRK2) and one gene related to the inflammation and cytokine system (P2RX7). We replicated association of the functional variant Val66Met of BDNF with BP and better performance in retention. The variants at the 5’ end of SLC6A4 also showed some evidence of association among males (p = 0.004), but the widely studied functional variants did not yield any significant results. A protective four-variant haplotype on P2RX7 showed evidence of association with BP and executive functions: semantic and phonemic fluency (p = 0.006 and p = 0.0003, respectively). Third, we analyzed 23 bipolar families originating from the North-Eastern region of Finland. A genome-wide scan was performed using the 6K single nucleotide polymorphism (SNP) array. We identified susceptibility loci at chromosomes 7q31 with a LOD score of 3.20 and at 9p13.1 with a LOD score of 4.02. We followed up both linkage findings in the complete set of 179 Finnish bipolar families. The finding on chromosome 9p13 was supported (maximum LOD score of 3.02), but the susceptibility gene itself remains unclarified. In the fourth part of the thesis, we wanted to test the role of the allelic variants that have associated with bipolar disorder in recent genome-wide association studies (GWAS). We could confirm findings for the DFNB31, SORCS2, SCL39A3, and DGKH genes. The best signal in this study comes from DFNB31, which remained significant after multiple testing corrections. Two variants of SORCS2 were allelic replications and presented the same signal as the haplotype analysis. However, no association was detected with the PALB2 gene, which was the most significantly associated region in the previous GWAS. Our results indicate that BP is heterogeneous and its genetic background may accordingly vary in different populations. In order to fully understand the allelic heterogeneity that underlies common diseases such as BP, complete genome sequencing for many individuals with and without the disease is required. Identification of the specific risk variants will help us better understand the pathophysiology underlying BP and will lead to the development of treatments with specific biochemical targets. In addition, it will further facilitate the identification of environmental factors that alter risk, which will potentially provide improved occupational, social and psychological advice for individuals with high risk of BP.

The brain serotonin transporter binding in young adults; methodological considerations and association with Bulimia Nervosa and acquired obesity

The neurotransmitter serotonin (5-HT) modulates many functions important for life, e.g., appetite and body temperature, and controls development of the neural system. Disturbed 5-HT function has been implicated in mood, anxiety and eating disorders. The serotonin transporter (SERT) controls the amount of effective 5-HT by removing it from the extracellular space. Radionuclide imaging methods single photon emission tomography (SPET) and positron emission tomography (PET) enable studies on the brain SERTs. This thesis concentrated on both methodological and clinical aspects of the brain SERT imaging using SPET. The first study compared the repeatability of automated and manual methods for definition of volumes of interest (VOIs) in SERT images. The second study investigated within-subject seasonal variation of SERT binding in healthy young adults in two brain regions, the midbrain and thalamus. The third study investigated the association of the midbrain and thalamic SERT binding with Bulimia Nervosa (BN) in female twins. The fourth study investigated the association of the midbrain and hypothalamic/thalamic SERT binding and body mass index (BMI) in monozygotic (MZ) twin pairs. Two radioligands for SERT imaging were used: [123I]ADAM (studies I-III) and [123I]nor-beta-CIT (study IV). Study subjects included young adult MZ and dizygotic (DZ) twins screened from the FinnTwin16 twin cohort (studies I-IV) and healthy young adult men recruited for study II. The first study validated the use of an automated brain template in the analyses of [123I]ADAM images and proved automated VOI definition more reproducible than manual VOI definition. The second study found no systematic within-subject variation in SERT binding between scans done in summer and winter in either of the investigated brain regions. The third study found similar SERT binding between BN women (including purging and non-purging probands), their unaffected female co-twins and other healthy women in both brain regions; in post hoc analyses, a subgroup of purging BN women had significantly higher SERT binding in the midbrain as compared to all healthy women. In the fourth study, MZ twin pairs were divided into twins with higher BMI and co-twins with lower BMI; twins with higher BMI were found to have higher SERT binding in the hypothalamus/thalamus than their leaner co-twins. Our results allow the following conclusions: 1) No systematic seasonal variation exists in the midbrain and thalamus between SERT binding in summer and winter. 2) In a population-based sample, BN does not associate with altered SERT status, but alterations are possible in purging BN women. 3) The higher SERT binding in MZ twins with higher BMIs as compared to their leaner co-twins suggests non-genetic association between acquired obesity and the brain 5-HT system, which may have implications on feeding behavior and satiety.

Studies on Causes of Multiple Sclerosis : From Genes to Transcriptome

Multiple sclerosis (MS) is the most common cause of neurological disability in young adults, affecting more than two million people worldwide. It manifests as a chronic inflammation in the central nervous system (CNS) and causes demyelination and neurodegeneration. Depending on the location of the demyelinated plaques and axonal loss, a variety of symptoms can be observed including deficits in vision, coordination, balance and movement. With a typical age of onset at 20-40 years, the social and economic impacts of MS on lives of the patients and their families are considerable. Unfortunately the current treatments are relatively inefficient and the development of more effective treatments has been impeded by our limited understanding of the causes and pathogenesis of MS. Risk of MS is higher in biological relatives of MS patients than in the general population. Twin and adoption studies have shown that familial clustering of MS is explained by shared genetic factors rather than by shared familial environment. While the involvement of the human leukocyte antigen (HLA) genes was first discovered four decades ago, additional genetic risk factors have only recently been identified through genome-wide association studies (GWAS). Current evidence suggests that MS is a highly polygenic disease with perhaps hundreds of common variants with relatively modest effects contributing to susceptibility. Despite extensive research, the majority of these risk factors still remain to be identified. In this thesis the aim was to identify novel genes and pathways involved in MS. Using genome-wide microarray technology, gene expression levels in peripheral blood mononuclear cells (PBMC) from 12 MS patients and 15 controls were profiled and more than 600 genes with altered expression in MS were identified. Three of five selected findings, DEFA1A3, LILRA4 and TNFRSF25, were successfully replicated in an independent sample. Increased expression of DEFA1A3 in MS is a particularly interesting observation, because its elevated levels have previously been reported also in several other autoimmune diseases. A systematic review of seven microarray studies was then performed leading to identification of 229 genes, in which either decreased or increased expression in MS had been reported in at least two studies. In general there was relatively little overlap across the experiments: 11 of the 229 genes had been reported in three studies and only HSPA1A in four studies. Nevertheless, these 229 genes were associated with several immunological pathways including interleukin pathways related to type 2 and type 17 helper T cells and regulatory T cells. However, whether these pathways are involved in causing MS or related to secondary processes activated after disease onset remains to be investigated. The 229 genes were also compared with loci identified in published MS GWASs. Single nucleotide polymorphisms (SNP) in 17 of the 229 loci had been reported to be associated with MS with P-value less than 0.0001 including variants in CXCR4 and SAPS2, which were the only loci where evidence for correlation between the associated variant and gene expression was found. The CXCR4 variant was further tested for association with MS in a large case-control sample and the previously reported suggestive association was replicated (P-value is 0.0004). Finally, common genetic variants in candidate genes, which had been selected on the basis of showing association with other autoimmune diseases (MYO9B) or showing differential expression in MS in our study (DEFA1A3, LILRA4 and TNFRSF25), were tested for association with MS, but no evidence of association was found. In conclusion, through a systematic review of genome-wide expression studies in MS we have identified several promising candidate genes and pathways for future studies. In addition, we have replicated a previously suggested association of a SNP variant upstream of CXCR4 with MS. Keywords: autoimmune disease, common variant, CXCR4, DEFA1A3, HSPA1A,gene expression, genetic association, GWAS, MS, multiple sclerosis, systematic review

Molecular basis of colorectal cancer predisposition

Colorectal cancer (CRC) is one of the most frequent malignancies in Western countries. Inherited factors have been suggested to be involved in 35% of CRCs. The hereditary CRC syndromes explain only ~6% of all CRCs, indicating that a large proportion of the inherited susceptibility is still unexplained. Much of the remaining genetic predisposition for CRC is probably due to undiscovered low-penetrance variations. This study was conducted to identify germline and somatic changes that contribute to CRC predisposition and tumorigenesis. MLH1 and MSH2, that underlie Hereditary non-polyposis colorectal cancer (HNPCC) are considered to be tumor suppressor genes; the first hit is inherited in the germline and somatic inactivation of the wild type allele is required for tumor initiation. In a recent study, frequent loss of the mutant allele in HNPCC tumors was detected and a new model, arguing against the two-hit hypothesis, was proposed for somatic HNPCC tumorigenesis. We tested this hypothesis by conducting LOH analysis on 25 colorectal HNPCC tumors with a known germline mutation in the MLH1 or MSH2 genes. LOH was detected in 56% of the tumors. All the losses targeted the wild type allele supporting the classical two-hit model for HNPCC tumorigenesis. The variants 3020insC, R702W and G908R in NOD2 predispose to Crohn s disease. Contribution of NOD2 to CRC predisposition has been examined in several case-control series, with conflicting results. We have previously shown that 3020insC does not predispose to CRC in Finnish CRC patients. To expand our previous study the variants R702W and G908R were genotyped in a population-based series of 1042 Finnish CRC patients and 508 healthy controls. Association analyses did not show significant evidence for association of the variants with CRC. Single nucleotide polymorphism (SNP) rs6983267 at chromosome 8q24 was the first CRC susceptibility variant identified through genome-wide association studies. To characterize the role of rs6983267 in CRC predisposition in the Finnish population, we genotyped the SNP in the case-control material of 1042 cases and 1012 controls and showed that G allele of rs6983267 is associated with the increased risk of CRC (OR 1.22; P=0.0018). Examination of allelic imbalance in the tumors heterozygous for rs6983267 revealed that copy number increase affected 22% of the tumors and interestingly, it favored the G allele. By utilizing a computer algorithm, Enhancer Element Locator (EEL), an evolutionary conserved regulatory motif containing rs6983267 was identified. The SNP affected the binding site of TCF4, a transcription factor that mediates Wnt signaling in cells, and has proven to be crucial in colorectal neoplasia. The preferential binding of TCF4 to the risk allele G was showed in vitro and in vivo. The element drove lacZ marker gene expression in mouse embryos in a pattern that is consistent with genes regulated by the Wnt signaling pathway. These results suggest that rs6983267 at 8q24 exerts its effect in CRC predisposition by regulating gene expression. The most obvious target gene for the enhancer element is MYC, residing ~335 kb downstream, however further studies are required to establish the transcriptional target(s) of the predicted enhancer element.

Insights into the molecular genetics of celiac disease : Applying family- and population-based strategies

Celiac disease, or gluten intolerance, is triggered by dietary glutens in genetically susceptible individuals and it affects approximately 1% of the Caucasian population. The best known genetic risk factors for celiac disease are HLA DQ2 and DQ8 heterodimers, which are necessary for the development of the disease. However, they alone are not sufficient for disease induction, other risk factors are required. This thesis investigated genetic factors for celiac disease, concentrating on susceptibility loci on chromosomes 5q31-q33, 19p13 and 2q12 previously reported in genome-wide linkage and association studies. In addition, a novel genotyping method for the detection of HLA DQ2 and DQ8 coding haplotypes was validated. This study was conducted using Finnish and Hungarian family materials, and Finnish, Hungarian and Italian case-control materials. Genetic linkage and association were analysed in these materials using candidate gene and fine-mapping approaches. The results confirmed linkage to celiac disease on the chromosomal regions 5q31-q33 and 19p13. Fine-mapping on chromosome 5q31-q33 revealed several modest associations in the region, and highlighted the need for further investigations to locate the causal risk variants. The MYO9B gene on chromosome 19p13 showed evidence for linkage and association particularly with dermatitis herpetiformis, the skin manifestation of celiac disease. This implies a potential difference in the genetic background of the intestinal and skin forms of the disease, although studies on larger samplesets are required. The IL18RAP locus on chromosome 2q12, shown to be associated with celiac disease in a previous genome-wide association study and a subsequent follow-up, showed association in the Hungarian population in this study. The expression of IL18RAP was further investigated in small intestinal tissue and in peripheral blood mononuclear cells. The results showed that IL18RAP is expressed in the relevant tissues. Two putative isoforms of IL18RAP were detected by Western blot analysis, and the results suggested that the ratios and total levels of these isoforms may contribute to the aetiology of celiac disease. A novel genotyping method for celiac disease-associated HLA haplotypes was also validated in this thesis. The method utilises single-nucleotide polymorphisms tagging these HLA haplotypes with high sensitivity and specificity. Our results suggest that this method is transferable between populations, and it is suitable for large-scale analysis. In conclusion, this doctorate study provides an insight into the roles of the 5q31-q33, MYO9B, IL18RAP and HLA loci in the susceptibility to celiac disease in the Finnish, Hungarian and Italian populations, highlighting the need for further studies at these genetic loci and examination of the function of the candidate genes.

Molecular features of colorectal cancer predisposition and progression

Both inherited genetic variations and somatically acquired mutations drive cancer development. The aim of this thesis was to gain insight into the molecular mechanisms underlying colorectal cancer (CRC) predisposition and tumor progression. Whereas one-third of CRC may develop in the context of hereditary predisposition, the known highly penetrant syndromes only explain a small fraction of all cases. Genome-wide association studies have shown that ten common single nucleotide polymorphisms (SNPs) modestly predispose to CRC. Our population-based sample series of around thousand CRC cases and healthy controls was genotyped for these SNPs. Tumors of heterozygous patients were analyzed for allelic imbalance, in an attempt to reveal the role of these SNPs in somatic tumor progression. The risk allele of rs6983267 at 8q24 was favored in the tumors significantly more often than the neutral allele, indicating that this germline variant is somatically selected for. No imbalance targeting the risk allele was observed in the remaining loci, suggesting that most of the low-penetrance CRC SNPs mainly play a role in the early stages of the neoplastic process. The ten SNPs were further analyzed in 788 CRC cases, 97 of which had a family history of CRC, to evaluate their combined contribution. A significant association appeared between the overall number of risk alleles and familial CRC and these ten SNPs seem to explain around 9% of the familial clustering of CRC. Finding more CRC susceptibility alleles may facilitate individualized risk prediction and cancer prevention in the future. Microsatellite instability (MSI), resulting from defective mismatch repair function, is a hallmark of Lynch syndrome and observed in a subset of all CRCs. Our aim was to identify microsatellite frameshift mutations that inactivate tumor suppressor genes in MSI CRCs. By sequencing microsatellite repeats of underexpressed genes we found six novel MSI target genes that were frequently mutated in 100 MSI CRCs: 51% in GLYR1, 47% in ABCC5, 43% in WDTC1, 33% in ROCK1, 30% in OR51E2, and 28% in TCEB3. Immunohistochemical staining of GLYR1 revealed defective protein expression in homozygously mutated tumors, providing further support for the loss of function hypothesis. Another mutation screening effort sought to identify MSI target genes with putative oncogenic functions. Microsatellites were similarly sequenced in genes that were overexpressed and, upon mutation, predicted to avoid nonsense-mediated mRNA decay. The mitotic checkpoint kinase TTK harbored protein-elongating mutations in 59% of MSI CRCs and the mutant protein was detected in heterozygous MSI CRC cells. No checkpoint dysregulation or defective protein localization was observable however, and the biological relevance of this mutation may hence be related to other mechanisms. In conclusion, these two large-scale and unbiased efforts identified frequently mutated genes that are likely to contribute to the development of this cancer type and may be utilized in developing diagnostic and therapeutic applications.

Genetics of T cell co-stimulatory receptors -CD28, CTLA4, ICOS and PDCD1 in immunity and transplantation

Co-stimulatory signals are essential for the activation of naïve T cells and productive immune response. Naïve T cells receive first, antigen-specific signal through T cell receptor. Co-stimulatory receptors provide the second signal which can be either activating or inhibitory. The balance between signals determines the outcome of an immune response. CD28 is crucial for T cell activation; whereas cytotoxic T lymphocyte associated antigen 4 (CTLA4) mediates critical inhibitory signal. Inducible co-stimulator (ICOS) augments cytokine expression and plays role in immunoglobulin class switching. Programmed cell death 1 (PDCD1) acts as negative regulator of T cell proliferation and cytokine responses. The co-stimulatory receptor pathways are potentially involved in self-tolerance and thus, they provide a promising therapeutic strategy for autoimmune diseases and transplantation. The genes encoding CD28, CTLA4 and ICOS are located adjacently in the chromosome region 2q33. The PDCD1 gene maps further, to the region 2q37. CTLA4 and PDCD1 are associated with the risk of a few autoimmune diseases. There is strong linkage disequilibrium (LD) on the 2q33 region; the whole gene of CD28 exists in its own LD block but CTLA4 and the 5' part of ICOS are within a same LD block. The 3' part of ICOS and PDCD1 are in their own separate LD blocks. Extended haplotypes covering the 2q33 region can be identified. This study focuses on immune related conditions like coeliac disease (CD) which is a chronic inflammatory disease with autoimmune features. Immunoglobulin A deficiency (IgAD) belongs to the group of primary antibody deficiencies characterised by reduced levels of immunoglobulins. IgAD co-occurs often with coeliac disease. Renal transplantation is needed in the end stage kidney diseases. Transplantation causes strong immune response which is tried to suppress with drugs. All these conditions are multifactorial with complex genetic background and multiple environmental factors affecting the outcome. We have screened ICOS for polymorphisms by sequencing the exon regions. We detected 11 new variants and determined their frequencies in Finnish population. We have measured linkage disequilibrium on the 2q33 region in Finnish as well as other European populations and observed conserved haplotypes. We analysed genetic association and linkage of the co-stimulatory receptor gene region aiming to study if it is a common risk locus for immune diseases. The 2q33 region was replicated to be linked to coeliac disease in Finnish population and CTLA4-ICOS haplotypes were found to be associated with CD and IgAD being the first non-HLA risk locus common for CD and immunodeficiencies. We also showed association between ICOS and the outcome of kidney transplantation. Our results suggest new evidence for CTLA4-ICOS gene region to be involved in susceptibility of coeliac disease. The earlier published contradictory association results can be explained by involvement of both CTLA4 and ICOS in disease susceptibility. The pattern of variants acting together rather than a single polymorphism may confer the disease risk. These genes may predispose also to immunodeficiencies as well as decreased graft survival and delayed graft function. Consequently, the present study indicates that like the well established HLA locus, the co-stimulatory receptor genes predispose to variety of immune disorders.

Nephrin-associated molecules in human glomerular diseases

The glomerular epithelial cells and their intercellular junctions, termed slit diaphragms, are essential components of the filtration barrier in the kidney glomerulus. Nephrin is a transmembrane adhesion protein of the slit diaphragm and a signalling molecule regulating podocyte physiology. In congenital nephrotic syndrome of the Finnish type, mutation of nephrin leads to disruption of the permeability barrier and leakage of plasma proteins into the urine. This doctoral thesis hypothesises that novel nephrin-associated molecules are involved in the function of the filtration barrier in health and disease. Bioinformatics tools were utilized to identify novel nephrin-like molecules in genomic databases, and their distribution in the kidney and other tissues was investigated. Filtrin, a novel nephrin homologue, is expressed in the glomerular podocytes and, according to immunoelectron microscopy, localizes at the slit diaphragm. Interestingly, the nephrin and filtrin genes, NPHS1 and KIRREL2, locate in a head-to-head orientation on chromosome 19q13.12. Another nephrin-like molecule, Nphs1as was cloned in mouse, however, no expression was detected in the kidney but instead in the brain and lymphoid tissue. Notably, Nphs1as is transcribed from the nephrin locus in an antisense orientation. The glomerular mRNA and protein levels of filtrin were measured in kidney biopsies of patients with proteinuric diseases, and marked reduction of filtrin mRNA levels was detected in the proteinuric samples as compared to controls. In addition, altered distribution of filtrin in injured glomeruli was observed, with the most prominent decrease of the expression in focal segmental glomerulosclerosis. The role of the slit diaphragm-associated genes for the development of diabetic nephropathy was investigated by analysing single nucleotide polymorphisms. The genes encoding filtrin, densin-180, NEPH1, podocin, and alpha-actinin-4 were analysed, and polymorphisms at the alpha-actinin-4 gene were associated with diabetic nephropathy in a gender-dependent manner. Filtrin is a novel podocyte-expressed protein with localization at the slit diaphragm, and the downregulation of filtrin seems to be characteristic for human proteinuric diseases. In the context of the crucial role of nephrin for the glomerular filter, filtrin appears to be a potential candidate molecule for proteinuria. Although not expressed in the kidney, the nephrin antisense Nphs1as may regulate the expression of nephrin in extrarenal tissues. The genetic association analysis suggested that the alpha-actinin-4 gene, encoding an actin-filament cross-linking protein of the podocytes, may contribute to susceptibility for diabetic nephropathy.

Characterization of LRRTM and NGR gene families: Expression and functions

Some leucine-rich repeat (LRR) -containing membrane proteins are known regulators of neuronal growth and synapse formation. In this work I characterize two gene families encoding neuronal LRR membrane proteins, namely the LRRTM (leucine-rich repeat, transmembrane neuronal) and NGR (Nogo-66 receptor) families. I studied LRRTM and NGR family member's mRNA tissue distribution by RT-PCR and by in situ hybridization. Subcellular localization of LRRTM1 protein was studied in neurons and in non-neuronal cells. I discovered that LRRTM and NGR family mRNAs are predominantly expressed in the nervous system, and that each gene possesses a specific expression pattern. I also established that LRRTM and NGR family mRNAs are expressed by neurons, and not by glial cells. Within neurons, LRRTM1 protein is not transported to the plasma membrane; rather it localizes to endoplasmic reticulum. Nogo-A (RTN4), MAG, and OMgp are myelin-associated proteins that bind to NgR1 to limit axonal regeneration after central nervous system injury. To better understand the functions of NgR2 and NgR3, and to explore the possible redundancy in the signaling of myelin inhibitors of neurite growth, I mapped the interactions between NgR family and the known and candidate NgR1 ligands. I identified high-affinity interactions between RTN2-66, RTN3-66 and NgR1. I also demonstrate that Rtn3 mRNA is expressed in the same glial cell population of mouse spinal cord white matter as Nogo-A mRNA, and thus it could have a role in myelin inhibition of axonal growth. To understand how NgR1 interacts with multiple structurally divergent ligands, I aimed first to map in more detail the nature of Nogo-A:NgR1 interactions, and then to systematically map the binding sites of multiple myelin ligands in NgR1 by using a library of NgR1 expression constructs encoding proteins with one or multiple surface residues mutated to alanine. My analysis of the Nogo-A:NgR1 -interactions revealed a novel interaction site between the proteins, suggesting a trivalent Nogo-A:NgR1-interaction. Our analysis also defined a central binding region on the concave side of NgR1's LRR domain that is required for the binding of all known ligands, and a surrounding region critical for binding MAG and OMgp. To better understand the biological role of LRRTMs, I generated Lrrtm1 and Lrrtm3 knock out mice. I show here that reporter genes expressed from the targeted loci can be used for maping the neuronal connections of Lrrtm1 and Lrrtm3 expressing neurons in finer detail. With regard to LRRTM1's role in humans, we found a strong association between a 70 kb-spanning haplotype in the proposed promoter region of LRRTM1 gene and two possibly related phenotypes: left-handedness and schizophrenia. Interestingly, the responsible haplotype was linked to phenotypic variability only when paternally inherited. In summary, I identified two families of neuronal receptor-like proteins, and mapped their expression and certain protein-protein interactions. The identification of a central binding region in NgR1 shared by multiple ligands may facilitate the design and development of small molecule therapeutics blocking binding of all NgR1 ligands. Additionally, the genetic association data suggests that allelic variation upstream of LRRTM1 may play a role in the development of left-right brain asymmetry in humans. Lrrtm1 and Lrrtm3 knock out mice developed as a part of this study will likely be useful for schizophrenia and Alzheimer s disease research.

Genetic studies on recurrent miscarriage

Recurrent miscarriage (RM) is defined as three consecutive pregnancy failures and is estimated to affect ~1% of couples trying to conceive. The cause of RM remains unknown in approximately 50% of cases. In this study, it was hypothesized that some of the underlying factors yet to be discovered are genetic. The aim was to search for mutations in genes AMN, EPCR, TM, and p53 known to cause miscarriage in mouse models and thereby find new genetic causes for unexplained miscarriages in humans. In addition, the mitochondrial genome was studied because mitochondria are involved in processes important in early development. Furthermore, sex chromosome characteristics suggested to underlie miscarriage were also studied. A total of 40 couples and 8 women with unexplained RM were collected for this study and screened for mutations in the candidate genes. Six interesting exonic or potential splice site disrupting variations were detected. However, their phenotypic effects cannot be determined without further investigations. Additionally, an association between the C11992A polymorphism of the p53 gene and RM was detected. The results indicate that women carrying the C/A or A/A genotype have a two-fold higher risk for RM than women with a C/C genotype. This strengthens the results of previous studies reporting that p53 sequence variations may cause miscarriage. The role of variation C11992A in embryonic development is, however, difficult to predict without further studies When screening the mitochondrial genome a heteroplasmic mtDNA variation was found in an unexpected high number of women, as heteroplasmic variations are reported to be rare. One novel variation and 18 previously reported polymorphisms were detected in the mitochondrial genome. Although the detected variations are likely to be neutral polymorphisms, a role in the aetiology of miscarriage cannot be excluded as some mtDNA variations may be pathogenic only when a threshold is reached. Recent publications have reported skewed X chromosome inactivation and Y chromosome microdeletions to be associated with RM. Therefore, these sex chromosome abnormalities in the context of RM were investigated. No associations between skewed X chromosome inactivation or Y chromosome microdeletions and RM in the Finnish patients were detected. Data on ancestral birthplaces of the patients were collected to study any possible geographic clustering, which would indicate a common predisposing factor. The results showed clustering of the birthplaces in eastern Finland in a subset of patients. This suggests a possibility of an enriched susceptibility gene which may contribute to RM.